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Relativistic fluctuating hydrodynamics 

Esteban Calzetta* 

lAFE and Physics Department, UBA, Buenos Aires, Argentina 
We derive the formulae of fluctuating hydrodynamics appropiate to a 

relativistically consistent divergence type theory, obtaining Landau - Lifshitz 

fluctuating hydrodynamics as a limiting case. 

I. INTRODUCTION 



Boltzmann equation [0], and then the hydrodynamic limit is taken by conventional means 



^ \ By applying the fluctuation - dissipation theorem [||] to the Navier - Stokes equation, 

\ Landau and Lifshitz developed long ago the standard theory of fluctuating hydrodynamics 

■ 01 HI- This theory is both appealing and successful, and it is consistent with the fluctuating 

OO \ hydrodynamics which is derived if fluctuation - dissipation theory is first applied to the 

O 

OO 

' 0- It has also been applied to black hole fluctuations and vacuum decay [^. The simplest 

' relativistic generalization of Navier - Stokes theory being the Eckart theory of dissipative 

^ ' fluids H (by simplicity, in this paper we shall not discuss the related Landau - Lifhsitz theory 
0], nor other 'first order' theories |^), it is only natural to provide a relativistic version of 

•rH . 

^ ' fluctuating hydrodynamics by applying fluctuation - dissipation theory to it; this step has 



been taken, and the results applied to such fields as cosmology fT^- Trouble is, Eckart 



theory has strong drawbacks as a relativistic theory, among them lack of stable solutions 
and acausal propagation of perturbations [|ri|]. This paper asks what changes must be 
introduced in Landau - Lifshitz fluctuating hydrodynamics if the Eckart theory is replaced 
by a suitable relativistic theory. 

Unfortunately, there is not a single causal version of relativistic real fluid dynamics 



as compelling as the Eckart and Landau - Lifshitz theories |]T2[ [|T3|, the main contenders 
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being the extended thermodynamics theories (ETTs) [|T4| |T5| and the divergence type 



theories (DTTs) |T8[, both of which draw some support from relativistic kinetic theory 



g pg pTf (see PI for ETTs, pi for DTTs). In particular the ETT framework has been 



extensively applied to cosmology . The relationship between these two approaches is not 
yet well understood. For concreteness, we shall adopt the DTT framework, which is more 
adapted to a rigorous statement of results. 

It ought to be said that the issues raised in this paper are to some extent academic, 
since there are strong reasons to believe that in any case the Eckart type theories are a 



good phenomenological description of any relativistic fluid ||2g] However, it is to 
be expected that small differences will become important when the theory is pushed to 
its limits, such as the relativistic theory of dissipative superfluids (see for a 



ETT approach to this problem, for a DTT one, and for a different perspective). 



the use of relativistic hydrodynamics as a phenomenological approximation to semiclassical 
cosmological models , or the emergence of hydrodynamic descriptions as a preferred set 



of decoherent histories |Q ||^. It is with these applications in mind that the following 
considerations were developed. 

The next section provides a review of relativistic fluid dynamics; its aim, of course, is 
not to substitute the classic introductions such as refs. ||,|^,|^, but rather to have all the 
relevant results in one place and notation. We then give a primer on Landau - Lifshitz 
fluctuation theory, and develop the subject of fluctuations in DTTs. 

II. RELATIVISTIC FLUID DYNAMICS 

A. Basic thermodynamics 

Consider a system described by extensive quantities like entropy, energy, volume, particle 
number and momentum S, U, V, N, P and intensive ones like temperature, pressure, chemical 
potential and velocity T,p,fi,v. 
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We build a covariant theory by adopting the following rules : 

a) Intensive quantities (T, p, fi) are associated to scalars, which represent the value of the 
quantity at a given event, as measured by an observer at rest with respect to the fluid. 

b) Extensive quantities {S, V, N) are associated to vector currents 5", u", N"- (we assume 
MTW conventions, c = 1, signature -+++, and latin indexes go from to 3 such that 
given a time like surface element dT^a = UadT,, then — X'^c/Sq is the amount of quantity X 
within the volume dH as measured by an observer with velocity n"'. If further the quantity X 
is conserved, then X.^ = 0. The quantity associated to volume is the fluid four-velocity, 
and obeys the additional constraint u"^ = —1. 

c) Energy and momentum are combined into a single extensive quantity and associated 
to a tensor current T"*. The energy current, properly speaking, is U"' = —T°'^Ub. 

The entropy current S"" is given by 

which we rewrite as 

where we introduced the affinity a = fi/T, the thermodynamic potential = p/?", and 
the inverse temperature vector (3°' = T'^u"'. The Gibbs - Duhem relation becomes 

TdS" = -diT^'^Ub) + pdu" - fidN" 

We now introduce the concept of a 'perfect fluid', as a system whose energy - momentum 
tensor takes the form T"'' = pu"'u''+pA"''', where A"* = u"'u^+g"'^ (observe that p must be the 
energy density seen by an observer moving with the fluid). Since we must have Uadu"' = 0, 
we conclude pdu"" — T"'^duf, = 0. Thus for a perfect fluid, the Gibbs - Duhem relation reads 

dS^ = -adN" - (3bdT^^ 

and 
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d^" = N^da + T'^'^dPh 
which is how the thermodynamic potential got its name. Equivalently 



The symmetry of the energy momentum tensor imphes that the thermodynamic potential 
is itself a gradient 

where Xp is the so called generating function. In a covariant theory, Xp ^lay only depend 
on the scalars a and PaP"", working out the corresponding derivatives, we conclude that 
necessarily A^" = nw", where n is the particle number density seen by a comoving observer. 

We conclude with a word on equilibrium states. Suppose the fluid departs from equi- 
librium by a fluctuation ^A^"", (5T"^, consistent with the conservation laws but otherwise 
arbitrary. Then the change in entropy production is 

But for a true equilibrium state the entropy must be stationary, and so we must have 
C(,a = P{b;a) = |12| (here and henceforth, brackets stand for symmetrization) . Thus the 
affinity must be constant, and the inverse temperature vector must be Killing. 

B. Eckart's theory of real fluids 

We wish now to describe a weakly dissipative fluid. Following Eckart, we shall base our 
description on the same set of variables than for an ideal fluid, namely, T"-^ and A^*^, only 
now the energy - momentum tensor is given as 

T"'' = pu^'u^ + pM"^ + t"^ (3) 

where it is assumed that the viscous stress t""^ obeys 



and p is defined to be the same function of p and n than in the ideal case. Assuming 
that the entropy current and thermodynamic potential are also the same, we get the entropy 
production rate 



SI 



(4) 



Now expand 



Pa;b = [-TbUa + TUa;b] 



ab I „ab 



(5) 



where 



to get 



SI = 



The principle of positive entropy production is satisfied if 



(6) 



(7) 



(8) 



(9) 



representing heat flux, bulk and shear viscosity, respectively. With these constitutive re- 
lations, energy - momentum conservation yields an straightforward, covariant generalization 
of the Navier - Stokes theory. 
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can be written directly in terms of the covariant derivatives of the inverse temperature 
vector as 



where 



nb riabcd n 

^ = -B P(c;d) 



(10) 



B 



abed 



AkT^P^'""^ + AriTP^"^" + QCbvTPs 



abed 



yabed 



(11) 



pabed j_ ^ ab ^ cd 



(12) 



ryabed 



^a^e^bd _^ ^a^d^be _^ ^b^e^ad _^ ^^b^d^ae 



(13) 



-yabed 
TT 



^ae ^bd _j_ ^ad ^ 



ad A be 



'_^ab ^ed 



(14) 



Observe that B 



abed 



5 , and that the P's are actually projection operators, that is, 



= P in all three cases, while the product of different P's vanishes. 

Eckart's theory is such a compelling generalization of non relativistic dissipative hydro- 
dynamics that it is a pity it doesn't work. The resulting equations allow for non causal 



propagation, and ipso facto all their solutions are unstable ||11 



While it is fairly clear that Eckart's theory (and the closely related Landau's theory as 
well) must be rejected, it is not clear at all what should be their replacement. The so - 
called Israel - Stewart or 'second order' type theories fl^ [|l^ ||16| perform much better 
with regards to both causality and stability, while keeping much of the appeal of the Eckart 
framework, but still lack a rigorous proof of consistency. 

In the following, we shall adopt instead the Geroch - Lindblom 'divergence type' de- 
scription of a relativistic real fluid [|I^ . The resulting theory is further removed from direct 
thermodynamic intuition than the Eckart proposal, but does allow for a rigorous proof of 
causality and stability. 
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C. Divergence type real fluids 



The failure of the Eckart approach to real fluids may be attributed to two unwarranted 
assumptions, namely, that the real fluid could be described within the same set of variables 
and with the same entropy current than its perfect counterpart. As a matter of fact, all 
that equilibrium thermodjTiamics suggests is that, whatever extra variables are brought in to 
describe the non equilibrium state, they must vanish in equilibrium, and the entropy current 
must match its equilibrium value up to first order in the deviations from equilibrium. 

According to Geroch and Lindblom |T^, description of a nonequilibrium state requires, 
besides the particle current and energy momentum tensor, a new third order tensor Aabc, 
obeying an equation of motion of divergence type 

Al;a = he (15) 

and 

Aabc = Aacb] ^fe =0; h = Q 

The entropy current is enlarged to read 

S^ = <^''- pf^T"'' - aN" - A^''Xbc (16) 

(ab vanishes identically in equilibrium, it is symmetric, and = 0. We further require 
entropy and the thermodynamic potential to be algebraic functions of their arguments. 

The condition that the conservation laws and Eq. ([T5|) imply positive entropy production 
demands that Eqs. (|l]) hold; in particular, the thermodynamic potential derives from a 
generating function as in Eq. (^. The thermodynamic potential is also allowed to depend 
on the new tensor (ab] actually 

^ = A^'^ (17) 

Thus the entropy production 
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St 



(18) 



We obtain positivity by making the former linearly dependent on the latter. 



The Eckart theory is actually a particular case of Geroch and Lindblom's [|I^ . Write the 
generating functional as 



XE = Xp + [^] CabU^U 



where Xp has the same form than for a perfect fluid. 
Recalling 



(19) 



Wb 



1 \ab . 



dT 

Wa 



'2„,a 



(20) 



we get the thermodynamic potential 



(21) 



The particle current is the same than for a perfect fluid, and the energy - momentum 
tensor is again of the form T"'^ = T!^^ + r"'', where the dissipative tensor 



ab /^abcd/- 



(22) 



C 



abed 



rp2 ^ pabcd _j_ pabcd _j_ _^ ^ab 



1 
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(in writing this equation, we have used that Ca = 0)- Observe that C is not symmetric. 
This equations allows us to define the viscous stresses in terms of Cafe? namely 

'A- 



„ab rp2 
71 = I 



\ac\bd _|_ \ad\bc _ gfc ^ cd 

~ 3 



cd 



(23) 



The new tensor is given by 



j^abc 



(24) 



and its divergence 



+ 



(25) 



T 



Observe that we can write 



A ecd /^Tcdab n 
^E;e = ^ P{a;b) 

where (J'^'^'^'^^ ^abcd 

Let us write Eqs. (|T5|, ^ and |2^) in shorthand as 



(26) 



T = CC 



These three equations ought to be equivalent to Eq. (|T0[) 



r = -BV(3 



To obtain this, we must provide a hnear relationship 



jab j~^abcd^ 



cd 



(27) 



or in shorthand 



/ = -DC 



(28) 



where 



D = C^B-^C 



(29) 



The inverse must exist, since B is positive definite. This equation may be inverted, to 
yield 
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B = CD-^C^ 



(30) 



where we assume that all matrices have inverses, even the non symmetric ones. 
Written in full 

4T3 



D 



abed 



BV 



T2 

Tjabcd _ pabcd , _ pabcd 



V 



(31) 



where 



^abcd 
ST 



(32) 



and the other projectors are defined in eqs. ([T^ and |T^). Observe that D is a symmetric 
operator. 

This concludes the setting of Eckart's theory in a DTT framework. 



D. Causality and divergence type theories 

While the Eckart theory belongs to the DTT class, it is not a 'good' theory, as we shall 
presently see. To investigate what conditions must a suitable DTT theory satisfy, we must 
discuss further the issues of causality and stability. 

The main advantage of the dissipative type theory framework is that the condition of 
causality may be expressed in a remarkably simple form. Let us introduce the symbol to 
denote the triad (a,/?",C''), A% the triad {N'',T''\ A"'"'), and Ib the triad (0,0,4b). Then 
the theory is summed up in the equations (cfr. Eqs. (P jr7| , p!8|JT5D ) 

4a _ 



The equations of motion can also be written as 
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where 



= M^cB = (33) 



Then the causahty condition is that the quadratic form M'^^Wa be negative definite for 
all future directed timelike vectors w"", or, equivalently, that for any displacement 6^^ from 
an equilibrium state, the vector = M'^q6(^^6('" be timelike and future oriented . 

In order to see the meaning of this condition, it is interesting to observe that under any 
such displacement, the change in the entropy current is precisely (55"° = —2(5", and therefore 
the change in entropy density, as seen by an observer with velocity Va, is 6s = 2vaQ"'. Thus 
the causality condition estates that the entropy will be reduced by any displacement from 
equilibrium. This goes beyond mere thermodynamic stability, which only requires entropy 
reduction at constant energy and particle number densities. 

For imperfect fiuids, it is necessary to consider the full quadratic form, that IS, as a 
functional of Sa, 5(3°" and ^C"*- K is easy to see that the Eckart theory cannot possibly meet 



the test, since a whole diagonal block is missing. Geroch and Lindblom [0 have suggested 
a simple way of constructing causal theories close to Eckart 's. The idea is to write down a 
generating functional of the form (cfr. Eq. (|19])) 



x = xi. + 7(^)g(w^c') (34) 

where Q is a positive definite scalar quadratic form on the (^"^.Then the thermodynamic 
potential 

= $^ + u^QT^^ + A"^T7|^ (35) 
di ou° 



and, provided 



d'-y 

5r» 



> 0, (36) 
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the second term dominates and guarantees the positivity of the corresponding vector Q". 
In what follows we shall assume one such extension of Eckart's theory has been adopted. 

This closes our review of causal relativistic hydrodynamics. We now proceed to our main 
subject, namely, the form the fluctuation - dissipation theorem takes in this context. 



III. FLUCTUATION THEORY 

A. The classical fluctuation - dissipation theorem 

The simplest possible application of fluctuation - dissipation theory relates to an homo- 
geneous system described by time - dependent variables Let 

Then the entropy production rate is given by 

S = -Xix' 

Following common usage, we refer to the X's as thermodynamic forces, and the x's 
as thermodynamic fluxes. Then the principle of positive entropy production on average is 
satisfied by posing a linear relationship 

x\t) = - fdt' j'^{t, t')Xj{t') + C{t) (37) 



where 7 is positive definite (Onsager's reciprocity principle further asserts that 7*-' = ±7-'* 
according to whether behaves as under time reversal or not, and we also assume 
causality). Under equilibrium, we have the average 

which follows from Einstein's formula relating entropy to fluctuations, and averages are 
time traslation invariant. Therefore 

x^x') + (x'x^) = 
12 



If the noise ^ is Gaussian, then 



x^c) = J dt' {ce{t') 



Prom the equations of motion, we conclude that for a given initial condition x^{t'), the 
solution is written in terms of a single Green function as 

x\t) = G% t')x={t') + dr G% r)e(r) 

where = 5]. Therefore 

5x'{t) _ 5x\t) 
WW) ~ ~Mf) 

(which expresses in symbols Onsager's insight that the regression of microscopic fluctu- 
ations follows the same rules than macroscopic ones) 
Thus 

(ee) = (^) [7^^' + 7^'] (39) 

which is our basic result. 

The relationship of forces to fluxes may be inverted to yield 

Xi = GijX^ + (40) 

Then it is immediate that 

{e,{t)e,{t')) = (\)[G,, + G,,] (41) 



.2 

We may also wish to introduce new fluxes y — Mx, conjugated to forces Y — M~^X. 
The relationship between fluxes and forces becomes 

y^RY + a; M-fM : a ^ 

Then 

(aa) = R 

This generalized fluctuation dissipation theorem will be relevant in what follows. 
We shall now apply this basic framework to Eckart's theory and to a causal relativistic 
theory. 
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B. Fluctuations in Eckart's theory 



In order to derive the spontaneous fluctuations of a Eckart type real fluid, we first 
regard a fluid theory as an instance of the above example, with the index i now running 
over a continuous as well as a discrete range. We also go back to the entropy production 
rate Eq. (§) and identify r"'' as a 'flux' multiplying the corresponding 'force' P{a-b) (which, 
as it should, vanishes in equilibrium) (since we shall presently derive these results from a 
DTT viewpoint, we only give here the somewhat sketchy original Landau argument p|; the 
derivation is spelled out in full in Fox and Uhlembeck 0]). We then postulate a Langevin - 
type relationship 



where 5^^\x,x') is the covariant four dimensional delta function. The deterministic part 
is chosen to reproduce the Eckart constitutive relationships. 

Since the coefficients are already symmetric, the fluctuation dissipation theorem suggests 



On the other hand, if we decompose the stochastic energy - momentum tensor as in Eq. 




(42) 




(43) 



(§) 



ab 



(44) 



we find 



(qX) = 2«:T2 A-5(^) {x, x') 
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(g:(7r,A^'^ + vrf)) = (45) 

Which generahze the Landau - Lifshitz formulae to the relativistic regime [ |10|1 . 
We are finally prepared to consider the issue of fiuctuations in a relativistically consistent, 
divergence type theory. 

C. Fluctuations in divergence type theories 

In a divergence type theory, the entropy fiux takes the form (cfr. Eq. (|l^)) 

The definition 

implies the entropy production 

na /-B A a 

'->;a — S ^B;a 

The Langevin - type equations of motion read 

A%,a = Ib+ Jb (46) 

Where Jb is the stochastic noise. So far, all formulae in this section have been ultralocal, 
meaning that and Ib are algebraic functions of the field variables at the same point (as 
opposed, e. g., to depending on their derivatives). 

A satisfactory theory must predict vanishing mean entropy production in equilibrium, 
namely 

(C^(x)/^(x)) + (c^(a;)jB(x)) = (47) 

However, because the coincidence limit of the correlation functions in the left hand side 
may not be well defined, we must impose this condition in a smeared form, or else, coupling 
it to an elementary causality consideration, request the stronger condition 

15 



C^(x)Jb(x')) + (C''(x)jb(x')) = (48) 

for every spacelike pair {x,x'). 

In the linear approximation, the first term yields 

(C^(x)/B(a;')) = /(B,c) (c''(a;)C^(x')) (49) 

We now appeal to the general theory above to make the assumption that the noise is 
Gaussian, and, since the dynamics is purely local, white, meaning that 

{jB{x)jc{x')) = ascS'^'\x,x') (50) 

where ctbc is a given matrix. Under these hypothesis we obtain 

dC^ix) 



C^(x)jb(x')) = CTBC- 



djcix') 

Our goal is to derive the matrix asc, and our result shall be that indeed ctbc = —I(b,c)- 
To arrive to this end, we must show that, at equal times, 

<^''^'^'<^'') - IS 

Let us multiply both sides of this equation by the non singular matrix 

Mab = ria- 



Where is the unit normal to some Cauchy surface containing both x and x' (henceforth, 
the surface). In the linear approximation 

MabC = ^ 

where = —UaS"- is the entropy density in an adapted coordinate system. In equilib- 
rium, we may apply Einstein's formula (cfr. Eq. (pHf)) above, to conclude 



Mab {(''{xX'^ix')) = -6^6{x, x') (52) 
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where S{x,x') is the three dimensional covariant delta function on the Cauchy surface. 
On the other hand 

djcix') djc{x') 

where A^{x) = —UaA'^^x). Let us write the equations of motion Eq. (|46|) as 

dA^^ix) , . , , 
-i^ + L = Ux) 

where L involves the field variables on the surface, but not their normal derivatives, and 
d/dt = n'^d/dx"-. Indeed 

'li^-m^y) (54) 

(at equal times). From Eqs. (|52|), (|53D and (0), we obtain the desired result 



{jaJb) = { 



5C^ 6C^ 



(55) 



Eqs. (|46|) and ( ^5]) are the basic equations incorporating fluctuations to the DTT frame- 



work, and are our main result. 

While our derivation does not apply to Eckart's theory, since the matrix M above fails 
to be non singular, we may consider it as a limiting case of causal DTTs and still apply the 
final result. Since = Ipa = 0, and in equilibrium /^^^ = Ic,i3 ~ 0' S^^ reassuring 
result that particle number and energy momentum conservation are not violated. For the 
remaining equation we have, in the shorthand notation of the previous section 

= VA = I + j = -DC + J 

in) = D 

The first equation is an algebraic equation for the variables Its solution reads 

C = Cold + Cs, where (old corresponds to the usual expression in terms of the gradients in 
temperature and velocity (which may themselves be stochastic) and Cs = D^^j is a purely 
stochastic component with self correlation D^^. For the random viscous stresses we find 
(cfr. Eq. (H)) 
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s = CC 



and therefore 

(ss) = CD-^C^ 

But according to Eq. (PP]), CD~^C^ = B, and these are just the Landau - Lifshitz 
fluctuation formulae Eq. (0). 

This deceptively simple derivation should not disguise the fact that the Eckart theory 
is pathological. The fluctuation formulae Eqs. (^) and (^) can only be expected to yield 
sensible results in the context of truly causal theories, like those described in Eqs. 



and pGl) . 



D. Fluctuations in causal theories 

We finally arrive at our designated goal, namely, the characterization of hydrodynamic 
fluctuations in a truly consistent relativistic dissipative theory. The problem is that we do 
not have a unique generalization of Eckart 's theory to rely on, but rather a wide family of 
seemingly plausible extensions. We shall therefore appeal to Occam's razor, and concentrate 
on the simplest possible example, namely, the family of generalizations of Eckart's theory 
introduced in equations (|3^ , |35| and ^6]) above. While we shall appeal to rather drastic 
approximations to keep things simple, it should be clear that this serves illustration purposes 
only, and does not affect the validity of the noise correlation formula Eq. (|55D above. 

A DTT theory is defined by its generating functional and the driving forces in the 
equations of motion Eq. (pIS]). Since the shortcomings of Eckart's theory concern only 
the first, it is simplest to assume that the driving forces of the full theory remain the same, 
namely, they are still given by Eqs. and |31]) above. The generating functional is modified 



by adding a new term depending on a positive definite scalar Q (cfr. eq. (§^)) and a function 
7 of temperature. The simplest choice for the latter is a scale - free power law, 7 ~ T". Since 
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the restriction eq. (^) excludes the exponent n = 1, we shall settle for the next available 
alternative, n = 2. 

Concerning Q, we notice that the Eckart's theory already contains such an scalar, namely 
the entropy production rate —I"'^(ab, and it is simplest to assume that both are related. We 
therefore seek a Q scalar with the structure (cfr. eq. (PT|), and note that this Q depends 
explicitly on T) 



Q 



^^'''^ {(PstO - ^ {(PvO + — {(PttO 



(56) 



_3Cbv f^T 1] 

The physical meaning of the new r coefficients will become clear soon; for the time being 
we only remark that all three are assumed to be positive. As a matter of fact, using the 



Eckart's theory results as a guide, we know that in the end, contrasting eqs. (p2|) and (pSj) 
with eqs. (^, || and||), we shall obtain PstC ~ (bv: Pv( ~ f^T and PttC ~ V- Now in actual 
applications, shear viscosity effects are often more important than either bulk viscosity or 
heat conduction (cfr. ref. Q]). We could therefore retain only the third term, although that 
would mean our Q would be merely non negative, rather than positive definite. 

We shall appeal to the "2^1 approximation" (see Appendix) to retain only the first of 
the two new terms in the thermodynamic potential, eq. (^), obtaining 

$» = $^ + TV I2Q + j (57) 



where the Eckart potential is still given by eq. (0). Assuming for simplicity that the 
r's are independent of a (we shall of course assume that they are simply constant) the 
new addition to the thermodynamic potential does not affect the expression for the particle 
current in terms of field variables. The energy momentum tensor gets a new term, quadratic 
in the C'^'s, but we know that this term is small and can be neglected, at least in a first run 
(see below). The real change comes in the non equilibrium current A"-^'^, which now reads 

4:Tbv 



{Psrcr - 7^ (Pvcr + — {Pttcy 



where A'^'^ is the Eckart result eq. (|2^) . Its divergence becomes 
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(58) 



A'. 



abc 



A abc I rp3 



4r. 



BV 



be 



+ R 



be 



(59) 



where (de = u°'Cde;a- The Eckart term is given in eq. (^), and it is independent of 
the C's- The remainder i?'"^ is a comphcated expression involving products of (^"'"s and 
derivatives of the temperature and four velocity; it is therefore of second order in departure 
from equilibrium and can be neglected within present accuracy. 

As a consequence of Eq. (|59D , the third set of equations of motion is no longer merely 
algebraic; it now reads (cfr. eqs. and |3lD ) 

4rBy 



A% + T 

^3 



{PstC) 



K 



{PsrCr - 



2kT 



V 



4 / „ ^\bc / j-i /■\bc I 1 / n /■\bc 



BV 



, rj. iPvcr + - {PttCY 



+ J 



be 



(60) 



Using the orthogonality properties of the P's it is possible to decouple these equations; 
the simplest case is when all r's are equal, whereby we simply get 



=1) 



BV jjabed 
ST 



abed 
TT 



Eedie Jcd 



} 



(61) 



We may also use eqs. and to transform eqs. ( [5D| ) into a set of Maxwell - Cattaneo 
equations for the viscous stresses In any case, the r's represent the characteristic times 
in which the (^"^'s (and therefore the viscous stresses as well) relax to their Eckart values. 
The apparition of new coefficients with the physical meaning of relaxation times is also 
characteristic of the ETT approach (cfr. |[T^|1^) and to this extent we may conjecture that 
both approaches are physically equivalent within this accuracy. 

We can see now how the theory works. Because the /"^^ are the same, the noise statistics 
has not changed, and in particular we still have Gaussian white noise. The energy momen- 
tum tensor has not changed either, at least as regards the expression of the viscous stresses 
in terms of the (ab^ (there is a systematic contribution coming from the nonvanishing ex- 
pectation value of the new quadratic terms which ought to be included, though). What has 
changed is that, if the noise - noise correlation is ultralocal, then the stress - stress correla- 
tion cannot be, because the solution to eq. (|60| ) will be necessarily non local in the noises. 
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Of course, we expect these correlations will decay exponentially on a time of order r, and in 
the approximation in which these characteristic times are neglected, we recover the Landau 
- Lifshitz results above. 
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V. APPENDIX 



In this appendix, we shall take a closer view into the approximations leading to Eq. 
(^). The starting point is the generating functional (cfr. Eq. 



X = XE + 7iT)QiT,u-X 



,a /-ab 



(62) 



where Q is a positive definite scalar quadratic form on the Then the thermodynamic 



potential (cfr. Eq. (|35D ) 



dT 



(63) 



If 7 = and Q is given by Eq. 



we get 



u 



^ {(PstO - ^ iCPvO + — {(PttO 



+ A-T7g (64) 



Recall also Eqs. (|32| , |I3| and |T^). 

The basic idea is that the timelike second term will dominate the spacelike third term, 
thus ensuring negativity of M'^^Wa, where Wa is an arbitrary future oriented timelike vector, 
and M^^ is the Hessian Eq. (|33|). If this is true, then the third term in Eq. ( |6^ may be 
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discarded. Since the second term is larger than the third roughly by a factor 2, we have 
called this the "2^1 approximation" . 

(This argument does not hold for the heat conduction terms, which in fact are linear in 
T; for them, we must rely directly on the argument below, or else assume that thc has a 
hidden temperature dependence.) 

Of course, it is not hard to compute the third term explicitly, but is is not very illu- 
minating either. Our real interest lies in finding the tensor A"'^'^ and its divergence; when 
we do this, we find that the terms proportional to u*^ contribute time derivatives of as 
well as several extra terms, called R""^ in the text, Eq. (|59|), while the third term will only 
contribute space derivatives and cross terms containing both the tensor (^"^ and derivatives 
of the temperature and fluid velocity. Thus, in physical terms, the approximation involved 
is that evolution in time is more important than inhomogeneity in space. In the end we 
recover essentially the ETT equations of motion, which have their own physical motivation, 
independent of the formal procedure above. 
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